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Abstract. We discuss the influence of tidal spin-orbit interactions on the 
orbital dynamics of close intermediate- mass X-ray binaries. In particular we 
consider here a process in which spin angular momentum of a contracting 
RLO donor star, in a synchronous orbit, is converted into orbital angular 
momentum and thus helps to stabilize the mass transfer by widening the 
orbit. Binaries which would otherwise suffer from dynamically unstable 
mass transfer (leading to the formation of a common envelope and spiral-in 
evolution) are thus shown to survive a phase of extreme mass transfer on a 
sub-thermal timescale. Furthermore, we discuss the orbital evolution prior 
to RLO in X-ray binaries with low- mass donors, caused by the competing 
effects of wind mass loss and tidal effects due to expansion of the (sub)giant. 



1. Introduction 

Tidal torques act to establish synchronization between the spin of the non- 
degenerate companion star and the orbital motion. Whenever the spin an- 
gular velocity of the donor is perturbed (by a magnetic stellar wind; or 
change in its moment of inertia due to either expansion or mass loss in 
response to RLO) the tidal spin-orbit coupling will result in a change in 
the orbital angular momentum leading to orbital shrinkage or expansion. 
We have performed detailed numerical calculations of the non-conservative 
evolution of ~ 200 close binary systems with 1.0 — 5.0 Mq donor stars and 
a 1.3 Mq accreting neutron star. Rather than using analytical expressions 
for simple polytropes, we calculated the thermal response of the donor star 
to mass loss, using an updated version of Eggleton's numerical computer 
code, in order to determine the stability and follow the evolution of the 
mass transfer. We refer to Tauris & Savonije (1999) for a more detailed 
description of the computer code and the binary interactions considered. 
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2. The orbital angular momentum balance equation 

Consider a circular^ binary with an (accreting) neutron star and a compan- 
ion (donor) star with mass Mns and M2, respectively. The orbital angular 
momentum is given by: Jorb = (-^ns -^2 1^) , where M = Mns + -^2 
and $7 = yjGMja^ is the orbital angular velocity. A simple logarithmic dif- 
ferentiation of this equation yields the rate of change in orbital separation: 

a qMns 0M2 MNS + M2 

a Jorb Mns M2 M ^ ' 

where the total change in orbital angular momentum can be expressed as: 

Jorh Jgwr ^ Jmh ^ Jls ^ Jml 

Jorb Jorb Jorb Jorb Jorb 

The first term on the right side of this equation governs the loss of Jorb due 
to gravitational wave radiation (Landau Sz Lifshitz 1958). The second term 
arises due to a combination a magnetic wind of the (low-mass) companion 
star and a tidal synchronization (locking) of the orbit. This mechanism of 
exchanging orbital into spin angular momentum is referred to as magnetic 
braking (see e.g. Verbunt Sz Zwaan 1981; Rappaport et al. 1983). 



2.1. TIDAL TORQUE AND DISSIPATION RATE 

The third term in eq.(2) was recently discussed by Tauris & Savonije (1999) 
and describes possible exchange of angular momentum between the orbit 
and the donor star due to its expansion or mass loss (note, we have neglected 
the tidal effects on the gas stream and the accretion disk). For both this 
term and the magnetic braking term we estimate whether or not the tidal 
torque is sufficiently strong to keep the donor star synchronized with the 
orbit. We estimate the tidal torque due to the interaction between the 
tidally induced flow and the convcctive motions in the stellar envelope by 
means of the simple mixing-length model for turbulent viscosity u = aHpVc, 
where the mixing-length parameter a is adopted to be 2 or 3, Hp is the local 
pressure scaleheight, and Vc the local characteristic convective velocity. The 
rate of tidal energy dissipation can be expressed as (Terquem et al. 1998): 




where the integration is over the convcctive envelope and is the orbital 
angular velocity, i.e. we neglect effects of stellar rotation. The radial and 

^This is a good approximation since tidal effects acting on the near RLO giant star 
will circularize the orbit on a short timescale of ~ 10'' yr, cf. Verbunt & Phinney (1995). 
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horizontal tidal displacements are approximated here by the values for the 
adiabatic equilibrium tide: 

. , 2 1 dir%) 

^' = f'P[d^) ^'=6-r^r- (4) 

where for the dominant quadrupole tide (l = m= 2) / = —GM2/{^a^). 
The locally dissipated tidal energy is taken into account as an extra energy 
source in the standard energy balance equation of the star, while the cor- 
responding tidal torque follows as: T = —{\/Q){dE/dt). 
The thus calculated tidal angular momentum exchange dJ = Tdt between 
the donor star and the orbit during an evolutionary timestep dt is taken 
into account in the angular momentum balance of the system. If the so cal- 
culated angular momentum exchange is larger than the amount required to 
keep the donor star synchronous with the orbital motion of the compact star 
we adopt a smaller tidal angular momentum exchange (and corresponding 
tidal dissipation rate) that keeps the donor star exactly synchronous. 

2.2. SUPER-EDDINGTON ACCRETION AND ISOTROPIC RE-EMISSION 

The last term in eq. (2) is the most dominant contribution and is caused by 
loss of mass from the system (sec e.g. van den Heuvcl 1994; Sobcrman et 
al. 1997). We have adopted the "isotropic rc-cmission" model in which all 
of the matter flows over, in a conservative way, from the donor star to an 
accretion disk in the vicinity of the neutron star, and then a fraction, f3 of 
this material is ejected isotropically from the system with the specific orbital 
angular momentum of the neutron star. If the mass-transfer rate exceeds 
the Eddington accretion limit for the neutron star /3 > 0. In our calculations 
we assumed /3 = max[0, 1 - MEdd/-^2] and MEdd = 1-5 x IO'^Mq yr~^. 



3. Evolution neglecting spin-orbit couplings 

Assuming Jg^r = ^mb = -^is = and jml/<^orb = [iq^ M2/{M2 {I + q)) one 
obtains easily analytical solutions to eq.(l). In Fig. 1 we have plotted 

_ain(a)^ 0_ 3^-^ 

51n(g) ^g + l ^^9(1-/3) + 1 ^ ' 

as a function of the mass ratio q = M2/M^^. The sign of this quantity is 
important since it tells whether the orbit expands or contracts in response to 
mass transfer (note dq <G). We notice that the orbit always expands when 
q <1 and it always decreases when q > 1.28 [solving dln{a)/d\n{q) = for 
(3 = 1 yields q = {I + Vyf)/4. ~ 1.28]. If /? > the orbit can stih expand 
for 1< g < 1.28. Note, d\n{a)/dln{q) = 2/5 at q = 3/2 independent of p. 
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Figure 1. — 91n(a)/91ri(g) as a function of q for X-ray binaries. The different curves 
correspond to different constant values of j3 in steps of 0.1. Tidal effects wore not taken 
into account here. A cross is shown to highlight the case of g = 1 or 91n(a)/91n(g) = 0. 
The evolution during the mass-transfer phase follows these curves from right to left since 
M2 and q are decreasing with time (though fi need not be constant). 

4. Results including tidal spin-orbit couplings 

In Fig. 2 we have plotted the orbital evolution of an X-ray binary. The 
solid lines show the evolution including tidal spin-orbit interactions and 
the dashed lines show the calculations without these interactions. In all 
cases the orbit will always decrease initially as a result of the large initial 
mass ratio [q = 4.0/1.3 ~ 3.1). But when the tidal interactions are included 
the effect of pumping angular momentum into the orbit (at the expense of 
spin angular momentum) is clearly seen. The tidal locking of the orbit acts 
to convert spin angular momentum into orbital angular momentum causing 
the orbit to widen (or shrink less) in response to mass transfer/loss. The re- 
lated so-called Pratt &: Strittmatter (1976) mechanism has previously been 
discussed in the literature (e.g. Savonije 1978). Including spin-orbit interac- 
tions many binaries will survive an evolution which may otherwise end up 
in an unstable common envelope and spiral-in phase. An example of this is 
seen in Fig. 2 where the binary with initial Porb = 2.5 days (solid line) only 
survives as a result of the spin-orbit couplings. The dashed line terminating 
at M2 ~ 3.0 Mq indicates the onset of a run-away mass-transfer process 
(M2 > 10~^ Mq yr~^) and formation of a common envelope and possible 
collapse of the neutron star into a black hole. In fact, many of the systems 
with 2.0 < M2/MQ < 5.0 recently studied by Tauris, van den Heuvel &; 
Savonije (2000) would not have survived the extreme mass-transfer phase 
if the spin-orbit couplings had been neglected. 
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Figure 2. Evolution of orbital separation as a function of donor star mass during the 
RLO phase in a binary with M2 = 4.0 Mq (X=0.70, Z=0.02, a=2.0), Mns = 1.3 Mq and 
^orb ~ 8.0 and 2.5 days, top and bottom lines respectively. The lifetime of these X-ray 
binaries are only tx = 1-2 and 2.1 Myr, respectively. The solid evolutionary tracks were 
calculated including tidal interactions and the dashed lines without. See text for details. 

The location of the minimum orbital separations in Fig. 2 are marked by 

arrows in the case of Porb = 8.0 days. Since the mass-transfer rates in such 
an intermediate-mass X-ray binary are shown to be highly super-Eddington 
(Tauris, van den Heuvel & Savonije 2000) we have /? 1. Hence in the 
case of neglecting the tidal interactions (dashed line) we expect to find the 
minimum separation when q = 1.28 (cf. Section 3). Since the neutron star 
at this stage only has accreted ~ 10""^ Mq we find that the minimum orbital 
separation is reached when M2 = 1.28 x 1.30 = 1.66 Mq. Including tidal 
interactions (solid line) results in an earlier spiral-out in the evolution and 
the orbit is seen to widen when M2 < 1.92 {q « 1.48). 

4.1. LOW-MASS DONORS AND PRE-RLO ORBITAL EVOLUTION 

For low-mass (< 1.5 Mq) donor stars there are two important consequences 
of the spin-orbit interactions which result in a reduction of the orbital 
separation: magnetic braking and expansion of the (sub)giant companion 
star. In the latter case the conversion of orbital angular momentum into 
spin angular momentum is a caused by a reduced rotation rate of the donor. 
However, in evolved stars there is a significant wind mass loss (Reimers 
1975) which will cause the orbit to widen and hence there is a competition 
between this effect and the tidal spin-orbit interactions for determining the 
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Figure 3. The changes of donor mass, M2 (full lines) and orbital period, Porb (dashed 
lines), due to wind mass loss and tidal spin-orbit interaetions, from the ZAMS until the 
onset of the RLO as a function of the initial orbital period of a circular binary. 



orbital evolution prior to the RLO-phase. This is demonstrated in Fig. 3. 
We assumed M2wind = — 4x 10~^^ ryRw L R2/M2 Mq yr~^ where the mass, 
radius and luminosity are in solar units and ?7rw is the mass-loss parameter. 
It is seen that only for binaries with -P^b^^ > days will the wind mass 
loss be efficient enough to widen the orbit. For shorter periods the effects 
of the spin-orbit interactions dominate (caused by expansion of the donor) 
and loss of orbital angular momentum causes the orbit to shrink. This result 
is very important e.g. for population synthesis studies of the formation of 
millisecond pulsars, since Porb in some cases will decrease significantly prior 
to RLO. As an example a system with M2 = 1.0 M©, Mns = 1-3 Mq and 

PZAMS ^ 3 Q J^yg Qj^ly J^g^^g pRLO ^ I Q (Jg^yg ^ J^g Qj^gg^ ^ J^g J^J^ Q 
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